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Abstract 

The role of the isoscalar hyperon A in probing the density dependence of the nuclear 
symmetry energy is studied in multi-A hypernuclei, hyperon-rich matter and neutron stars 
in relativistic models. Relationships between the properties of three types of objects and 
the neutron thickness in 208 Pb are established with respect to the isoscalar-isovector cou- 
pling that modifies the density dependence of the symmetry energy. The exotic isotopes 
far from the neutron drip line can be stabilized by filling in considerable A hyperons. 
The difference of the binding energy of multi-A hypernuclei from different models is at- 
tributed to different symmetry energies. The isovector potential together with the neutron 
thickness in multi-A hypernuclei investigated is very sensitive to the isoscalar-isovector 
coupling. The large sensitivity of the A hyperon fraction to the isoscalar-isovector cou- 
pling occurs at about 2 ~ 3po in beta equilibrated hyperon-rich matter. In neutron stars 
with hypcronization, an on-off effect with respect to the isoscalar-isovector coupling exists 
for the neutron star radius. 

PACS: 21.60.-n, 21.80.+a, 26.60.+C, 27.60.+j 

Keywords: Hypernuclei, bulk matter, neutron stars, relativistic mean field model 



'Email: jiangwz02@hotmail.com, Also at Center of Theoretical Nuclear Physics, National Laboratory of 
Heavy Ion Accelerator, Lanzhou 730000, China 



1 



The symmetry energy is pervasively related to the descriptions of nuclear structures and 
reactions in laboratories and astrophysics [1, 2, 3]. However, the density dependence of the 
symmetric energy is still poorly known[4, 5] to date. Therefore, theoretical uncertainties exist 
in the understanding of a large number of properties of neutron stars and heavy nuclei[6, 7], 
nuclei far from the (3 stability [8, 9, 10], heavy ion collisions[5, 11] and consequently probable 
hadron-quark phase transition[12], etc. Many factors that influence the isovector potential 
directly or indirectly can modify the behavior of the density dependence of the symmetry 
energy. Recently the density dependence of the symmetry energy has been extensively ex- 
plored through the inclusion of the isoscalar-isovector coupling terms in relativistic mean 
field (RMF) models[6, 7, 8, 13]. These new terms enable one modify the neutron skin of 
heavy nuclei without changing the charge radius and a variety of ground-state properties 
that are well constrained experimentally. On the other hand, the charge density distribu- 
tions for some light exotic nuclei near drip lines that are not determined yet are sensitive to 
the isoscalar-isovector couplings[10, 14]. 

Hypernuclei that are produced through strangeness transfer reactions mostly contain 
only one hyperon, and just a few double A hypernuclear events were observed[15, 16, 17]. 
Through heavy-ion collisions it may possibly create an environment which is favorable to 
the formation of metastable exotic multihypernuclear objects (MEMOs) such as multi-A 
hypernuclei. Schaffner et.al. showed that multi-A hypernuclei were more strongly bound 
than normal nuclei using the RMF model[18, 19]. Recently, Jevorsek et.al. set limits on 
the existence and properties of exotic objects such as strangelets and MEMOs, and provided 
limits as low as 1CT 7 for M/Z up to 120 through a reanalysis of related data[20]. In bulk 
matter, strangeness can be formed by virtue of the strong interaction, and hyperons may be 
important constituents of neutron stars. Many important effects of hyperonization in neutron 
stars have been found in the past (for a review,see[21]). For instance, hyperonization reduces 
the maximum mass of neutron stars as much as 3/4M [22, 23]. Schaffner et.al. found that 
a phase transition to hyperon matter may even possibly occur in neutron stars [24]. In this 
work, it is the aim to constrain the density dependence of the symmetry energy through 
the investigation of properties of multi-A hypernuclei, hyperon-rich matter and neutron stars 
with hyperonization. We will perform the investigation in RMF models considering the 
isoscalar-isovector coupling terms that soften the symmetry energy at large densities. 
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The modification in MEMOs and hyperon-rich matter due to the isoscalar-isovector cou- 
plings is usually considered to be small due to the relatively weak meson-hyperon coupling. 
As pointed out in Refs.[18, 19] that the MEMOs may reach a very high baryon density 
(2.5 ~ 3po) like strangelets. In hyperon-rich matter, the appearance of A hyperons occurs at 
about twice normal density. Since the higher density relates to the larger isoscalar potentials 
(or fields), the isoscalar-isovector couplings will induce large modifications. In hyperon-rich 
matter, we will focus on the relation between the A fraction and the isoscalar-isovector cou- 
pling. As the determination of properties of neutron stars with hyperonization needs the 
equation of state of hyperon-rich matter, we may search for the observable consequence re- 
lating to the density dependence of the symmetry energy. As examples of MEMOs, we will 
consider multi-A hypernuclei based on the Ca isotope far off the normal neutron drip line at 
the fixed isospin. In this way, we can observe the extension of the drip line, and explore the 
density dependence of the symmetry energy by shifting the isoscalar field. 

In this study, additional strange mesons a* (i.e. fo, 975MeV) and <fi (1020MeV) are also 
included in RMF models to describe the strong AA attraction as in Refs.[19, 24, 25]. The 
effective Lagrangian density is given as follows 

£ = t/j B [i-f^ - M B + g aB a - g^sl^ - OpBl^nb^ 

+ ^k a ^ dU ^ - 4 (1 + t ^ a ^ b - u ^ ^ & o ) 

+\(d^a - m 2 a a 2 ) - \f, v F^ + 

-\B, V B^ + \mlb,X ~ ^W" + £y (1) 
and Cy is for the strange meson-hyperon interactions and free fields of strange mesons 
Cy = ^y{9**yo* -g^Yl^W + \(d^*d»a* -ml«a* 2 ) 

4 W " VWiW" ~ d M + (2) 

where ipB, c, to, and bo are the fields of the baryon, scalar, vector, and charge-neutral isovector- 
vector mesons, with their masses MB,m ff ,m u , and m p , respectively. is the field of the 
photon. gis(i = <t,uj,p) and f wB are the corresponding meson-baryon couplings. T3 is the 
third component of isospin Pauli matrix for nucleons and T3 = for the A hyperon. r c is a 
constant relating to the baryon charge. F^, B^, and are the strength tensors of the oj, 
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p mesons and the photon, respectively 

Fpv = d p uj„ - d u uj p , B pv = d p b 0l/ - d v b 0fl , A pv = d p A v - d v A p (3) 
The self-interacting terms of a, w-mesons and the isoscalar-isovector ones are in the following 

-4g 2 pN (A s g 2 aN a 2 + Kgl N u^)b X (4) 

For the hyperon-rich matter, the chemical equilibration is established on weak interactions of 
baryons and leptons as in Ref.[22, 26]. We study chemically equilibrated and charge neutral 
hyperon-rich matter including baryons (N, A, E) and leptons (e, p) without E hyperons, since 
E hyperons do not appear[24, 25] according to the fact that the isoscalar potential changes 
sign in the nuclear interior and becomes repulsive based on recent E _ atomic data [27]. For 
the finite system, we consider the spherical case. The solution detail can be referred to 
Refs.[28, 29], and it is not repeated here. The pairing correlation of nucleons in nuclei is 
considered in the BCS approximation as in Refs.[14]. 

We perform calculations with the NL3 [30] and S271[6] parameter sets. These two models 
share the same binding energy per nucleon (16.24MeV) and incompressibility (271MeV) at 
the same Fermi momentum /cf = 1.30fm _1 (pQ = 0.148fm~ 3 ). The symmetry energy is given 
by[7] 

k 2 a 2 k 3 

° 4 " 6E* F + 3n 2 mf (5j 



with Ep = \]k 2 F + M^ 2 and m* 2 = m 2 +Sg 2 N (g 2 aN K s a 2 +g 2 lN k v uol). The density dependence 
of the symmetry energy is different between two models because of different g p N, fn p and Mjy- 
(M* N = 0.59M N for the NL3 and M* N = 0.7M N for the S271). 

The lo and p meson couplings to the hyperons are given by the SU(3) relations: g UJ \ = 
2/3g w at = 2^=, <7 p a = and g p -= = g p N- The a meson coupling to the hyperons is determined 
by reasonable hyperon potentials in nuclear matter: = — 30MeV and = — 20MeV. 

The coupling constants g^ and g^ are taken to satisfy the SU(6) relation: g^s/g^N = 
2g<j)A/gu>N = — 2\/2/3, and the g a *A is fitted to improve the AA interaction matrix element 
A5aa[19]. In Ref.[24], the g a *A is taken in a range close to the g a N- We take the ratio 
of the scalar coupling constant to be g^h/gaN = 0.76 in NL3 with a potential for the A 
hyperon in S matter = — 55.6MeV. In S271, the same potential is taken to obtain the 



ratio g a *k/9(jN = 0.77. In this way, the A1?aa of ^Be is improved from 0.38 MeV to 2.2 
MeV in NL3, and from 0.2 MeV to 3.9 MeV in S271, much close to the extracted value 4 ~ 5 
MeV. The g a *s is determined by the S hyperon potential in E matter = — 40MeV[19]. 
The ujNN tensor coupling is vanishing, and the wAA tensor coupling is small but adjusted 
to simulate the vanishing spin-orbit splitting for the A hyperon observed in \ 6 0. 

For simplicity, we perform calculations for various isoscalar-isovector coupling A v 's, and 
the A s is set as zero. For a given coupling A v , we follow Ref.[6, 8] to readjust the pNN 
coupling constant g p N so as to keep an average symmetry energy fixed as 25.7 at hp = 1.15 
fm _1 . In doing so, it was found in Ref.[6] that the binding energy of 208 Pb is nearly unchanged 
for various A v 's. Below, we will see how the total binding energy of multi-A hypernuclei is 
modified by the A v . 



■ S271 




Figure 1: Binding energy per baryon for various Ca+A's systems as a function of the 
strangeness fraction f s = ua/ub at A v = 0. The symbols from left correspond to the 
hyperon number 0, 2, 8, 20, 28, 40, 50, 64, and 82, respectively. Except for results marked 
by the triangles, the strange mesons are considered. The solid and dashed curves in the inset 
are for the symmetry energy in NL3 and S271, respectively. 

The binding energy per baryon for various 102 Ca+A's systems is displayed in Fig.l. 
Metastable minima appears with the rise of the A hyperon number n\. The additional AA 
attraction, simulated by the strange mesons, produces much larger binding energy at larger 
f s and shifts the minima outwards. The difference of binding energies of multi-A hypernuclei 
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in NL3 and in S271 enlarges explicitly with the rise of n\. Since both models share the 
same saturation properties and take the same potential for the A hyperon, this enlargement 
is attributed to the different density dependence of the symmetry energy With the rise of 
riA, the central density of multi-A hypernuclei increases. For instance, the central density 
in 102 Ca+50A is about 1.43/?o in NL3 which is much larger than 1.06/?o in 102 Ca. The more 
different symmetry energy at the larger density domain, shown in the inset in Fig.l, thus 
accounts for the large difference of the binding energy at larger f s . 




Figure 2: Two-neutron separation energy for 104 Ca+A systems vs. f s . 

Fig.2 shows that the two-neutron separation energy of 104 Ca+A's systems increases from 
the negative to the positive with the rise of n\. As the n\ increases, the nucleon potential is 
enhanced by the attraction between hyperons and nucleons, being able to fill in more neutrons 
under the continuum. The neutron drip line can thus be largely extended. For instance, the 
drip line of Ca isotope can be extended from the neutron number 50 [31] to 82 as long as 
about 50 A's are filled in. The dependence of the separation energy on the isoscalar-isovector 
coupling A v reduces approximately with the increase of n\. In Refs.[18, 19], a large class 
of MEMOs are constructed by adding multi-hyperons in stable nuclei, while here we aim 
at the extension of the neutron drip line by filling in A hyperons in unstable isotopes. As 
multi-A hypernuclei may exert their roles in extending the neutron drip line and binding 
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baryons more strongly than in normal nuclei, they will potentially influence the astrophysical 
nucleosynthesis and the abundances of elements. To search for multi-A hypernuclei with 
the riA up to 50 is possible but certainly difficult [20]. Here, we just devote ourselves to 
investigating the sensitivity of properties of such exotic systems to the density dependence 
of the symmetry energy. 

The sensitivity of the total binding energy to the A v is different for various 102 Ca+A's 
systems. For 102 Ca+A's systems that are unstable against the neutron emission at small tja 
and unstable against the A hyperon emission at very large n\, the total binding energies can 
be modified by as much as tens of MeV for various A v 's. For stable 102 Ca+A's systems whose 
binding energies are situated on or very close to the minima, shown in Fig.l, the modification 
to the total binding energies by various A v 's is quite small. For instance, the total binding 
energy for 102 Ca+50A in NL3 is just modified by up to 2MeV with various A v values. 

A data to data correlation of the neutron thicknesses between 208 Pb and stable multi-A 
nuclei can be established as in Refs.[8, 10]. The neutron thickness of these stable multi-A 
nuclei is very sensitive to the A v . For instance, the uncertainty of the neutron thickness for 
102 Ca+40A and 102 Ca+50A is about 3 times that for 208 Pb. The proton density distribution 
in these exotic system is also sensitive to the A v and one third of the uncertainty of the 
neutron thickness comes from the protons. The large uncertainty of the neutron thickness is 
attributed to the addition of the A hyperons, which will be elucidated below. 

Usually, one modifies the isospin dependent potential considering the shift of the isospin 
and the isovector field directly. Here, we modify it by enhancing the isoscalar field (u field) 
through filling in A's without changing the isospin of the system, which is fulfilled by con- 
sidering the isoscalar-isovector coupling. In Fig.3, the isovector potentials are displayed for 
102 Ca and 102 Ca+50A in NL3. In S271, similar results for these two systems can be obtained. 
The isovector potentials for various A v 's in the exotic isotope are largely shifted by filling in 
considerable A hyperons. The uncertainty of the isovector potential against the A v in 102 Ca 
is much larger in the larger radius region, compatible with the understanding that 102 Ca is 
totally diffusive and unstable. Compared to results for 102 Ca, the sensitivity of the isovector 
potential to the A v in 102 Ca+50A is shifted to the central region, and at the large radius the 
convergence occurs, which is favored by the extension of the neutron drip line. The results 
shown in Fig.3 for two exotic systems are related to different density profiles. The central 
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Figure 3: The isovector potential for Ca and Ca+50A in NL3 vs. radius. The curves 
from above (near the arrow) are obtained with various A v 's (equably from 0.04 to 0.00) which 
give the neutron thicknesses (in fm) in 208 Pb as denoted in ascending order. 

density in 102 Ca+50A is about one third higher than that in 102 Ca. The deviation of the 
symmetry energy for various A v 's enlarges with the rise of the density beyond the normal 
density or even much lower density[4, 6, 7]. At the central region in 102 Ca+50A, the isovector 
potential is more sensitive to the A v than that in 102 Ca. The isoscalar potential can also be 
shifted to a certain degree by the A v , whereas for the normal nuclei the shift is negligible. 
The sensitivity of the isovector potential to the A v is responsible for the uncertainty of the 
neutron thickness. 

The A hyperon constituent is important in astrophysical bulk matter to explore the effects 
of the density dependence of the symmetry energy. We plot the A and H~ fractions defined by 
fy = Py/pb in Fig. 4 with the S271 based on a self-consistent calculation. The large sensitivity 
of the hyperon fractions to the A v occurs at about 2 ~ 3po where the A hyperon starts to 
appear. Since the density as high as 2 ~ 3po is accessible through the heavy ion collisions 
with the colliding energy about lGeV per nucleon, the detection related to the hyperons, 
especially the A's, in heavy ion collisions is helpful to probe the density dependence of the 
nuclear symmetry energy. The starting density for the A appearance is model dependent. For 
instance, it is about 1.8po m NL3 and 2.2/>o in S271. As shown in the inset in Fig.4, a shift 
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Figure 4: Hyperon fractions in S271 vs. baryon density in hyperon-rich matter. The curves 
from above are for different values of A v (equably from 0.00 to 0.05) that give the neutron 
thicknesses in 208 Pb as denoted in descending order. The inset shows the isovector potential 
with hyperons and without. 

of the isovector potential starts to appear with the A appearance, and the isovector potential 
differs more from that without hyperonization with the increase of the density. In multi-A 
nuclei far from the normal drip line, the isovector potential is shifted by filling in A's under 
the fixed isospin, while here the shift comes from the change of the isospin and the isoscalar 
field together. The convergent extent with respect to the A v exists for the hyperon fraction 
in the whole density domain. This indicates that it is favorable to include the isoscalar- 
isovector coupling terms for the model itself. At very high densities, the convergence of 
hyperon fractions occurs. This is attributed to the fact that the isospin is largely suppressed 
with the drop of the neutron fraction at high densities. The isovector potential is thus small 
compared to the isoscalar potential, so that the isoscalar potential is little changed through 
the isoscalar-isovector coupling. This implies that the hyperstars[24] whose core is solely 
composed of hypermatter at high densities will be insensitive to the A v . 

The large reduction of the maximum mass of neutron stars by hyperonization, firstly 
pointed out by Glendenning (see Ref.[21] and references therein), can be reproduced here. In 
NL3, hyperonization in the core reduces the maximum neutron star mass by about 0.9M Q 
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Figure 5: Radii of neutron stars as a function of the neutron thickness in 08 Pb. 



from about 2.8M , and by O.75M from about 2.4M in S271. Compared to models in which 
hyperons are absent, hyperonization in the neutron star whose mass is near the maximum 
mass increases the central baryon density by as much as lpo in the region r < 5km. For a 
1.4 solar-mass neutron star, the role of hyperonization is quite limited because of relative low 
core densities. Using the Oppenheimer-Volkoff equations for bulk matter in beta equilibrium, 
we thus investigate neutron stars with relative large masses: 1.8M with the NL3 model and 
1.62M with the S271 model. In Fig. 5, the radii of neutron stars with respective masses 
in two models are displayed as a function of the neutron thickness for 208 Pb. The radius R 
of neutron star without hyperonization increases with the neutron thickness in 208 Pb for a 
given parameter set, which is the same as in Ref.[7]. For neutron stars with hyperonization, 
the trend is similar, whereas the correlation between the R and the neutron thickness in 
208 Pb almost vanishes in a large region of the neutron thickness in 208 Pb. This property is 
consistent with the convergent extent shown in Fig.4. It indicates that for a neutron star with 
hyperonization the radius R nearly just depends on the on-off switch rather than concrete 
values of the isoscalar-isovector coupling in a given model. An accurate measurement of the 
neutron radius for 208 Pb at the Jefferson Laboratory[32] determines the density dependence 
of the symmetry energy at low densities, as pointed out in Ref.[7]. The measurement of 
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the radius of the neutron star with hyperonization can discriminate two regions of the A v 's 
to which the R is either very sensitive or insensitive at high densities. Various information 
through separate measurements may be combined together to constrain the isoscalar-isovector 
coupling that modifies the density dependence of the symmetry energy in the whole density 
region. 

In summary, we have studied the roles of the isoscalar A hyperons in exploring the density 
dependence of the symmetry energy in multi-A hypernuclei, hyperon-rich matter and neutron 
stars with RMF models. Relationships between the properties of three types of objects and 
the neutron thickness in 208 Pb are established with respect to the isoscalar-isovector coupling 
that modifies the density dependence of the symmetry energy. The exotic isotopes far from 
the neutron drip line can be stabilized by filling in considerable A hyperons, and the neutron 
drip line is substantially extended outwards. The difference of the binding energy of multi-A 
hypernuclei from different models is attributed to different symmetry energies. The isovector 
potential together with the neutron thickness in multi-A hypernuclei situated on the new 
drip line that is far from the normal drip line is very sensitive to the isoscalar-isovector 
coupling. In beta equilibrated hyperon-rich matter, we have investigated the sensitivity of the 
hyperon fractions to the isoscalar-isovector coupling based on the self-consistent calculation 
for hyperons and nucleons. The large sensitivity of the A hyperon fraction to the isoscalar- 
isovector coupling occurs at about 2 ~ 3po- The properties of neutron stars are investigated. 
An on-off effect with respect to the isoscalar-isovector coupling exists for the radius of neutron 
stars with hyperonization. Besides the detection of A hyperons in the heavy ion collisions, 
the measurements of the neutron star radius and neutron thickness in 208 Pb could help set 
up the constraints on the density dependence of the symmetry energy. 

Acknowledgement 

This work is partially supported by the National Natural Sciences Foundation of China under 
grant No.10405031. The author thanks Prof. M. Di Toro for helpful discussions. 

References 

[1] H. A. Bethe, Rev. Mod. Phys. 62, 801 (1990) 



11 



[2] I. Bombaci, in Isospin Physics in Heavy-Ion Collisions at Intermediate Energies, edited 
by B.A.Li and W. Udo Schroder (Nova Science Publishers, Inc., New York, 2001), p35 

[3] A.W. Steiner, M. Prakash, J.M. Lattimer, P.J. Ellis, Phys. Rept. 411 (2005) 325 

[4] B.A.Brown, Phys. Rev. Lett. 85(2000)5296 

[5] B.A.Li, Phys. Rev. Lett. 85(2000)4221; ibid.88(2002)192701 

[6] C.J.Horowitz, J.Piekarewicz, Phys. Rev. Lett.86(2001)5647 

[7] C.J.Horowitz, J.Piekarewicz, Phys. Rev. C64(2001)062802(R); ibid,C66(2002)055803 
[8] B.G.Todd, J.Piekarewicz, Phys. Rev. C67(2003)044317 
[9] Y.B. Wei Y.G. Ma W.Q. Shen, et.al. Phys.Lett.B586(2004)225 
[10] W.Z.Jiang, Y.L.Zhao, Phys. Lett. B617(2005)33 

[11] L.W.Chen, V.Greco, C.M.Ko, B.A. Li, Phys. Rev. Lett. 90(2003)162701 

[12] M. Di Toro, A. Drago, T. Gaitanos, V. Greco, A. Lavagno, nucl-th/0602052 

[13] J.Carriere, C.J.Horowitz, J.Piekarewicz, Astrophys. J 593(2003)463 

[14] W.Z.Jiang, Y.L.Zhao, Z.Y.Zhu, S.F.Shen, Phys. Rev. C72 (2005)024313 

[15] S.Aoki, et al Prog. Theor. Phys. 85(1991)1287 

[16] J.K.Ahn, et al Phys. Rev. Lett. 87(2001)132504 

[17] J.K.Takahashi, et al Phys. Rev. Lett. 87(2001)212502 

[18] J.Schaffner, C. B.Dover, A. Gal, C.Greiner, H.Stdcker, Phys. Rev. Lett. 71(1993)1328 

[19] J.Schaffner, C.B.Dover, A. Gal, C.Greiner, D.J. Millener, and H.Stdcker, Ann. Phys. 
(NY.) 235 (1994) 35 

[20] D.Javorsek II, E.Fischbach and D.Elmore,Phys. Rev. D67(2003)034015 

[21] N.K.Glendenning, Phys.Rept.342(2001)393 

[22] N.K. Glendenning, Astrophys. J. 293 (1985) 470 

12 



[23] N.K. Glendenning, S.A. Moszkowski, Phys. Rev. Lett. 67 (1991) 2414 

[24] J.Schaffner-Bielich, M.Hanauske, H.Stocker, W.Greiner, Phys. Rev. Lett. 
89(2002)171101 

[25] J.Schaffner-Bielich, A. Gal, Phys. Rev. C62(2000)034311 

[26] M.Prakash, I.Bombaci, M.Prakasha, P.J. Ellis, J.M. Lattimer and R.Knorren, Phys. 
Rept. 280(1997)1 

[27] J.Mares,E.Friedman, A. Gal, and B.K.Jennings, Nucl. Phys. A594(1995)311 

[28] B.D.Serot, J.D.Walecka, Adv. Nucl. Phys. 16(1986)1; Int.J.Mod.Phys. E6(1997)515 

[29] P.Ring, Prog. Part. Nucl. Phys. 37(1996)193 

[30] G.A.Lalazissis, J.K6nig, P.Ring, Phys. Rev. C55(1997)540 

[31] J.Meng, H. Toki, J.Y.Zeng, S.Q.Zhang, and S.G.Zhou, Phys. Rev. C 65 (2002) 041302(R) 

[32] R.Michaels, P.A.Souder, G.M.Urciuoli, spokespersons, Jefferson Laboratory Experiment 
E-00-003 



13 



